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SUMMARY 


The  potential  environmental  impact  of  supersonic  operations 
includes  damage  to  conventional  and  unconvenrional  structures  bv 
sonic  boom  overpressures.  The  assessment;  of  such  damage 
generally  requires  dedicated  flyovers  for  each  site  at  great 
expense.  In  this  report  we  show  that,  structural-acoustic 
techniques  may  be  used  to  help  provide  such  assessments  in  their 
absence.  These  techniques  allow  for  both  airborne  and 
groundborne  propagation  paths  and  also  account  for  diffraction. 
Tests  are  described  whereby  transfer  functions  relating 
structural  response  to  sonic  boom  overpressure  may  be  obtained 
using  a  broad  band  stationary  acoustic  source  and  approoriate 
data  processing  to  substitute  for  the  boom.  Further, 
structural-acoustic  reciprocity  is  invoked  allowing  these 
transfer  functions  to  be  determined  by  measuring  the  sound 
radiated  by  the  structure  when  driven  mechanically. 

Additionally,  it  is  demonstrated  that  state-of-the-art  general 
purpose  finite  element-boundary  element  structural-acoustic 
computer  codes,  which  inherently  satisfy  reciprocity,  now  provide 
a  powerful  analytical  tool  for  predicting  these  transfer 
functions  for  complex  site  geometries  and  constructions. 

Finally,  it  is  suggested  that  (1)  empirical  validation  of  the 
above  techniques  be  pursued  by  comparison  with  data  taken  at  a 
series  of  sites/structures  with  actual  supersonic  flyovers,  and 
(2)  upon  successful  completion  of  (1),  the  necessary  procedures 
be  codified  to  provide  a  new  structural  damage  assessment  tool 
for  those  responsible  for  site  planning  and  development  and 
operations  planning  in  vulnerable  environments. 


Iv 
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INTRODUCTION 


The  potential  environmental  impact  of  the  U.S.  Air  Force  or 
commercial  supersonic  operations  extends  to  humans  and  other 
animals  and  structures,  both  conventional  and  unconventional. 
Damage  may  be  psychological  or  physiological  in  the  first  case 
and  cosmetic  or  structural  in  the  second.  The  source  of  such 
damage  is  the  pressure  field  of  the  sonic  "boom".  Performing  an 
environmental  damage  assessment  typically  requires  dedicated 
aircraft  flyovers  for  each  site,  at  great  expense.  This  study 
investigates  a  technique  for  providing  such  assessments  in  their 
absence,  and  focuses  on  damage  to  structures.  The  approach  uses 
structural-acoustic  tests  in  the  reciprocal  mode.  It  is  based  on 
the  assumption  that  the  boom  pressure  levels  are  in  the  linear 
range  in  the  immediate  vicinity  of  the  subject  structures.  Wirh 
this  approach  a  stationary  vibration  source,  an  array  of 
microphones  and  digital  post  processing  of  the  measured  signals 
substitute  for  actual  flyovers.  In  Section  2.1  below  we  briefly 
review  current  procedures  for  relating  damage  assessment  to 
measurements.  A  brief  background  of  structural-acoustic 
reciprocity  is  provided  in  Section  2.2.1  and  in  2.2.2  we  outline 
the  general  procedure  for  the  proposed  application. 

2 . 1  STRUCTURAL  DAMAGE  ASSESSMENT 

Sonic  boom  induced  damage  to  structures  is  generally  envisioned 
in  terms  of  a  maximum  stress.  That  is  damage  occurs  if  and  when 
a  supersonic  operation  causes  a  stress  somewhere  in  the  structure 
to  exceed  a  material  strength.  For  example,  plaster  walls  and 
windows  are  particularly  susceptible.^  However  stress  is  not  a 
readily  measured  quantity.  As  a  consequence  surrogate  measures 
coupled  with  some  level  of  analytical  modelling  is  required.  For 
example  the  theory  of  linear  elasticity  allows  one  to  compute 
stress  from  strain  measurements  using  rosette  strain  gages  and 
estimates  of  material  elastic  moduli.  Peak  stresses  may  also  be 
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(1) 


predicted  from  measurements  of  structural  surface  pressures  or 
motions  by  cascading  transfer  functions,  analytically  or 
empirically  determined.  For  example, ^ 

Opk  =  Pa(>?p*/Pa)  =  P//(Pa/p//)(f5s/p 
with  the  dynamic  amplification  factor 

DAf  =  a^,/3,  (2) 

where 

0 pi,  =  peak  dynamic  stress  at  a  given  location  x  on  the 
structure  in  response  to  Pa 

dj  ==  stress  at  x  in  response  to  statically  applied  pressure 
with  magnitude  p^. 

P ff  =  free  field  acoustic  pressure 
Pa  =  pressure  on  the  structure 

The  first  ratio  in  Eg.  1  accounts  for  the  effects  of  acoustic 
diffraction  around,  and  reflections  from,  nearby  obstacles 
including  the  structure  in  question.  In  the  absence  of  such 
effects  the  ratio  is  unity.  The  second  ratio  characterizes  the 
induced  stress  in  the  structure  by  the  boom  pressure,  ignoring 
dynamic  effects  which  are  accounted  for  separately  in  the  DAF. 

In  Eg.  1  it  is  generally  the  case  that  the  quantities  p,-/  and 

perhaps  (Pa/P//)  are  measured  during  a  flyover  while  the  other 
functions  are  estimated  analytically. 


An  alternative 

formulation^  defines 

(^pyPa)  =  (ap*/t'p*)(hp,/Pa) 

(3) 

and  therefore 

=  Pa(0p*/I''pO(^'p*/Pa) 

(4) 

2 


where 

1/  =  eeak  structural  velocity  at  x  in  response  to  the 

applied  pressure.  The  advantage  of  Eq.  4  is  that,  under  a 
reasonably  broad  range  of  circumstances  the  ratio  ( a  p*. /  K  p*. )  may  be 
deterained  from  impedance-like  structural  measurements.^'^ 
Therefore  becomes  a  function  of  a  directly  measurable 

quantity,  possibly  but  not  necessarily,  during  a  flyover.  For 
example,  it  may  be  shown  that  this  ratio  takes  on  the  simple  form 

where  Pp  is  the  surface  mass  density  and  Cp  the  compressional 
sound  speed  in  the  material  with  K=\  for  compressional  waves  in 
a  rod,  K  =  i3  for  a  rectangular  beam  in  flexure  and  1.12<A:<2.0 
for  flexing  plates  with  aspect  ratios  ranging  broadly  from  0 . 1  to 

1.0. 

To  this  point  the  applicable  damage  assessment  metric  has  been 
taken  to  be  peak  stress,  either  directly  or  indirectly  via  peak 
velocity.  However  in  some  cases  stress  is  not  necessarily  the 
most  useful  measure.  As  an  example  consider  the  rattling  of 
brie— a— brae  articles,  and  in  particular  the  onset  of  rattling 
which  may  provide  a  more  practical  criterion  than  actual 
breakage. 

Assume  that  the  rattling  of  bric-a-brac  articles  follows 
immediately  the  onset  of  their  sliding  relative  to  some  sonic 
boom  induced  base  motion.  Then,  assuming  a  horizontal  base 
acceleration,  cib,  and  a  one  dimensional  translational  model, 
slippage  is  avoided  provided  that  ah<[xg,  where  [i.  is  an 
appropriate  coefficient  of  static  friction  and  g  the  acceleration 
of  gravity.  Typically  p.  <  1 ,  for  example  0.4— 0.7  for  stone  on 
stone  and  0.25-0.5  for  wood  on  wood.  Therefore  this  yields  a 
more  stringent  criterion  than  had  we  considered  a  vertical  base 
motion  and  rattling  associated  with  loss  of  contact,  for  which 
the  criterion  is  ab<g.  To  account  for  tipping  one  must 
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intrcduca  the  center  of  gravity  and  the  moment  of  inertia  of  the 
article.  However,  here  too  the  criterion  may  be  expressed  in 
terms  of  the  base  acceleration. 

For  perspective,  a  sampling  of  empirically  based  damage  criteria 
for  the  peak  stress  and  structural  velocity  of  conventional  and 
unconventional  structures  is  shown  in  Tables  1  and  2  along  with  a 
qualitative  ranking  of  damage  suscepribility  by  structure  type 
for  historical  structures  in  Table  3.5  Interestingly,  the  levels 
shown  in  Tables  1  and  2  are  by  and  large  lower  than  those  given 
in  Ref.  6  for  "parricle  velocities  that  are  the  threshold  of 
possible  damage  for  both  commercial  structures  and 
residential-type  structures",  viz.,  4  and  1  in/sec  respectively. 

The  distinction  between  conventional  and  unconventional,  or 
historical,  or  as  suggested  in  Ref.  2  irreplaceable,  structures 
is  not  fundamental  to  the  validity  of  structural-acoustic 
testing.  Rather  the  issues  relate  to  practicality,  for  examole 
testing  accessibility,  and  the  precision  of  the  damage 
assessment,  for  example  owing  to  uncertainties  in  material 
properties  associated  with  aging.  The  latter  is  likely  in 
evidence  in  the  spread  of  data  shown  in  Tables  1  and  2. 

In  theory  structural-acoustic  reciprocity  is  compatible  with  all 
of  the  above  approaches  to  damage  assessment.  In  practice, 
however,  the  reciprocal  testing  of  complex  states  of 
stress/strain  is  not  deemed  feasible  and  we  will  confine 
ourselves  to  the  measurement  of  the  other  metrics,  viz.,  surface 
pressures  and  structural  accelerations,  velocities  and 
displacements,  to  be  used  in  conjunction  with  appropriate 
transfer  functions. 
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oung  S  Modulus  <6)  Esiimalcd  based  on  Peak  Strain  >  I6U  M  in/in  (Siagg,  cl  al..  19B4).  H  =  1.72  x  10°  psi  and  Ks  =  1.3^. 

Esiimaied  based  on  o  ■  K  (a/h)^  PS  K(a^)^  »  1 16  for  4  fi  x  8  fl  brick  wall  test  specimen. 

(8)  Yokel,  eiaJ..  1971,  esiimaied  stress  is  aciuall)'  computed  modulus  of  ivpiure  assuming  panial  end  fixiiy. 
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(')  Peak  velocity  of  slmctiire 

Peak  velocity  of  ground  at  base  of  structure 
As  cited  in  Siskind,  et  al.,  19SI)b.  Ap|)cndix  A 


Table  3.  Ranking*  of  Historical  Structures  According  to  Susceptibility 
to  Damage  from  Aircraft  Noise  (Table  4.1  of  Ref.  5) . 


Type  of  Structure 

Sonic  Boom 

Subsonic  Jet 

Heavy  Helicopter 

Historic  Sites 

Windows 

6 

2 

2 

Wood  frame,  plaster 

3 

3 

1 

Wood  frame,  wood  panels 

8 

7 

7 

Adobe 

7 

10 

8 

Masonry,  stone 

13 

13 

13 

Brick 

2 

9 

10 

Prehistoric  Sites 


Masonry/Stone  -  roof  intact 

4 

1 

3 

Adobe  -  roof  intact 

5 

4 

5 

Masonry/stone  -  no  roof 

11 

5 

6 

Adobe  -  no  roof 

10 

8 

9 

Seismically-sensitive  Areas 

Avalanche  -  loose  snow 

1 

6 

4 

Early  American  pictographs, 

12 

11 

11 

petroglyphs,  caves 

Avalanches  -  slab 

9 

12 

12 

Landslide  areas 

14 

14 

14 

*Rank  order  1  is  most  susceptible,  14  is  least  susceptible. 
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2.2  STRUCTURAL-ACOUSTIC  RECIPROCITY 


2.2.1  Background 

The  concept  of  reciprocity  for  linear  structural-acoustic  sysrems 
dates  back  to  the  mid-nineteenth  century.  Helmholtz's 
Reciprocity  Principle  (1360)  states  that  an  acoustic  signal 
remains  invariant  when  the  locations  of  a  compact  source  and 
receiver  are  interchanged  in  a  homogeneous  acoustic  medium 
containing  rigid  boundaries.  This  result  was  generalized  by 
Rayleigh'^  to  include  more  realistic  situations,  for  example, 
non-compact,  i.e.,  extended  sources  and  receivers,  the  presence 
of  elastic  scattering  boundaries  and  dissimilar,  coupled  acoustic 
media.  He  also  was  able  to  show  that  steady  flow  could  be 
accommodated  provided  its  direction  is  reversed  in  the  reciprocal 
situation,  a  result  further  explored  by  Lyamshev.®'^  On  the 
other  hand  it  has  been  shown  that  reciprocity  is  violated  by 
systems  encompassing  gyroscopic  forces  and  certain  classes  of 
dissipative  constraints, ^  factors  not  believed  to  be  crucial  to 
the  proposed  application. 

The  statement  of  reciprocity  in  the  field  of  mechanics  is  known 
as  Betti's  theorem.  Lyamshev  and  others  have  shown  that 
reciprocity  remains  valid  for  steady-state  structural-acoustic  or 
elasto-acoustic  systems  where  the  structure  is  driven  either 
mechanically  or  acoustically . ^ ^ The  dynamics  of  these  systems 
may  involve  all  forms  of  wave  motion,  e.g.,  compression,  shear, 
flexure,  torsion  etc.  Chertock  showed  that  reciprocity  was 
satisfied  for  the  nonsteady  or  transient  problem  as  well,  while 
investigating  the  response  of  a  submerged  complex  structure  to  an 
underwater  explosion. 

Structural-acoustic  reciprocity  has  been  applied,  extensively  and 
successfully,  to  a  wide  range  of  practical  problems.  These 
include  the  study  of  structures  submerged  underwater  where  the 
coupling  between  the  structure  and  the  acoustic  medium  is  strong. 
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and  structuras  vibrating  in  air  where  the  coupling  is  generally 
weak.  The  calibration  of  aiicrophones  and  hydrophones  is  perhaps 
the  earliest  application. Another  utilization  is  the 
evaluation  of  the  structure-borne  noise  component  of  propeller 
induced  aircraft  cabin  noise. Of  particular  note  is  the  work 
of  Ten  Wolde,  et  al  who  have  achieved  considerable  success  with 
their  reciprocity  experiments  aimed  at  diagnosing  parallel 
structural  propagation  paths  from  vibrating  ship-board  equipment 
to  sound  radiation  in  the  water. 

2.2.2  Proposed  Applicarion 

Our  generic  problem  is  pictured  in  Fig.  1.  An  aircraft 
travelling  supersonically  produces  a  sonic  "boom",  which  at  some 
distance  from  the  source  propagates  linearly  as  an  acoustic 
pulse.  The  pressure  field  impinges  on  the  ground  and  other 
obstacles,  i.e.,  structures,  at  an  angle  a  =  sin'^  ( 1 /.-W)  with 
M  =  U/c  where  U  is  the  flight  speed,  c  the  in  air  sound  speed  and 
M  the  effective  Mach  number.  In  the  vicinity  of  a  structure  the 
pressure  field  is  partially  diffracted  and  reflected.  In 
addition  it  may  transmit  energy  into  any  structure  either 
directly  or  indirectly  after  propagating  along  the  ground, 
generally  in  the  form  of  Rayleigh  surface  waves.  All  such 
effects  may  produce  stress  on  the  structure  and  therefore 
potentially  cause  damage.  And  all  these  effects  are  accounted 
for  with  our  reciprocal  approach  which  we  now  describe. 

We  refer  to  the  geometry  of  Fig.  2a  as  the  "direct"  problem  and 
for  convenience  it  will  be  analyzed  in  the  frequency  rather  than 
time  domain.  Specifically,  we  assume  that  a  distant  harmonic 
acoustic  source,  of  volume  velocity  is  responsible  for 

the  pressure  field  incident  on  the  structure.  Although  the 
pressure  field  near  the  source  may  be  highly  nonlinear  we  take 
the  field  in  the  vicinity  of  the  structure  to  be  linear.  This 
source  is  our  far  field  representation  of  the  sonic  boom 
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Acoustic  source 


Fig  2a.  Direct  Structural-Acoustic  Geometry 


Acoustic  receiver 


p(R;o) 


Fig  2b.  Reciprocal  Structural-Acoustic  Geometry 


generation.  It  is  defined  as  the  equivalent  acoustic  source  that 
reproduces  the  incident  pressure  field  on  the  structure  for  a 
given  flyover,  Pyy(i?':cu).  Assuming  a  simple  monopole  source 

Q(i?:co)  =  [4T[i/?'-/?|  /(-ia)p)]exp[-/A:  \R'-R\]Pij{R'\uo)  (6) 

where  p  is  the  air  density,  and  k  =  uo/c  the  acoustic  wavenumber  at 
frequency  co.  The  structure  responds  linearly,  but  otherwise 
generally  as  an  elastic  or  viscoelastic  fabrication,  with 
velocity  v(x;uo'). 

We  now  pose  the  "reciprocal"  situation  sketched  in  Fig.  2b.  A 
harmonic  force  f(x;co)  is  applied  to  our  structure  at  x.  The 
structure  vibrates  in  response,  producing  the  acoustic  field 
p(^R;(x>).  Invoking  reciprocity^^ 


v(x;(x>)  p(R;ai) 

Q(R;(X))  F(x:(ju) 

(7) 

v(^x  ;uj)  =  Q(R  ■,(x>')p(R  -.(jo)/  F(x’,uo) 

(8) 

If  we  now  interpret  Q(/?:co)  as  the  frequency  spectrum  of  a  time 
signature  then  the  time  series  of  the  response  becomes 


v(x  if)  =  (2ji')  ^  J 


(9) 


and  the  peak  velocity  may  be  determined.  In  other  words,  the 
peak  velocity  of  our  structure  at  any  and  all  locations  in 
response  to  an  incident  pressure  field  may  be  determined  by 
mechanically  driving  the  structure  at  that  location  and  measuring 
the  associated  radiated  pressure.  Wind  or  other  anomalous 
propagation  factors  are  of  little  consequence  since  range 
requirements  are  at  most  the  order  of  acoustic  wavelengths. 
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The  only  requirement  on  the  spectrum  of  the  mechanical  source  is 
adequate  signal  to  noise  over  the  frequency  band  of  interest.  It 
need  not  mimic  the  spectrum  of  the  boom  since  the  reciprocal 
relationship  being  invoked  in  Eq.  7  is  in  the  form  or  a  spectral 
ratio.  For  example  any  particular  boom  spectrum  may  be 
synthesized  by  post  processing  impulse  data  using  an  instrumenred 
hammer.  We  now  show  that  the  force  levels  required  for 
reciprocal  testing  are  moderate. 

The  pressure  radiated  to  range  R  hy  a.  point  driven,  effectively 
infinite,  thin  elastic  panel  may  be  estimated  from 

I  P,^^CR-,(x>)  i=  ip/p  JF/2jiRh 

where  F  is  the  spectral  level  of  the  force,  (p/pp)  the  ratio  of 
air  to  panel  mass  densities  and  h  is  the  panel  thickness.  For 
illustrative  purposes,  taking  a  1  lb.  (4.4N)  force  on  a  2.5  cm 
thick  wood  panel  of  specific  gravity  0.5  at  a  nominal  range  of 
Im,  we  obtain  the  (frequency  invariant)  pressure  spectral  level 
of  roughly  70  dB  re  0.0002  [ibar.  The  corresponding  drive 
acceleration  of  this  panel  is  given  by 

I  a(cu)  1=  (/3/4)ajf /ppCph^ 

where  Cp  is  the  sound  speed  in  the  panel  material.  Letting 
Cp  =  3. 5x  10^ vi/ sec,  this  force-yields  a  drive  acceleration  of  -60  + 
20  log  f  (Hz)  dB  re  Ig.  For  perspective,  the  acceleration  of  our 
effectively  infinite  panel  to  a  normally  incident  acoustic  wave 
of  magnitude  p,  is  given  by  the  mass  law, 

I  a  1=  2pi/pp/i 

and  for  instance  a  1  psf  (47.3  N/m^)  incident  pressure  yields  an 
acceleration  level  of  -2.5  dB  re  Ig. 
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3 .  STRUCTURAL-ACOUSTIC  PREDICTIVE  MODELS  FOR  THE  RESPONSE  OF 
STRUCTURES  TO  SONIC  BOOMS 

The  excitation  of  a  structure  by  the  pressure  field  of  a  sonic 
boom  may  be  both  airborne  and  seismic  and  may  result  from  direct 
and/or  diffracted  pressure  components.  In  the  following  three 
subsections  we  demonstrate  structural-acoustic  reciprocity  for 
each  of  these  phenomena,  treated  separately  for  the  sake  of 
simplicity.  Both  analytical  models  (Sections  3.1  and  3.3)  and  a 
numerical  model  (Section  3.2)  are  presented  to  display  the  broad 
spectrum  of  predictive  tools  that  are  now  available,  particularly 
with  general  purpose  structural-acoustic  finite  element-boundary 
element  (FE/BE)  computer  codes. 

3 . 1  THE  ILLUMINATED  AIRBORNE  PATH 

To  analyze  the  airborne  path  for  illuminated  structural  elements 
we  consider  the  idealized  geometry  shown  in  Fig.  3. a.  A  plane 
wave  is  incident  on  a  simply  supported  rectangular  thin  plate. 

The  plate,  which  may  be  orthotropic  to  simulate  the  effects  of 
framing,  responds  in  flexure.  Vibrational  energy  dissipation  is 
accounted  for  by  means  of  a  structural  loss  factor  (ri)  which  is 
incorporated  into  the  elastic  modulus  (F),  making  it  complex, 
i.e.,  F=^F(l-ir|).  In  computing  the  drive  pressure  on  the  plate 
it  is  assumed  to  be  baffled.  The  approach  is  to  first  obtain  the 
solution  in  the  frequency  domain  for  a  pure  impulse.  The 
solution  for  a  given  waveform,  or  pressure  time  history,  is  then 
obtained  by  performing  the  necessary  convolution,  or  inverse 
(Fourier)  transform.  Next,  it  is  shown  that  this  solution  may  be 
interpreted  reciprocally  as  the  sound  pressure  radiated  by  the 
plate  when  driven  mechanically.  Finally,  numerical  results  are 
presented  for  an  illustrative  example  and  related  to  peak  stress. 


Acoustic  source 


Fig  3a.  Simply  Supported  Plate  Insonified  By  Acoustic 
Source:  Direct  Geometry 


3.1.1  Plate  Response  To  Incident  Acoustic  Wave 


The  harmonic  displacement  of  an  orthotropic,  rectangular,  simoly 
supported  thin  plate  to  the  imposed  point  pressure  f'(uj)6(Xs)6(y ^.) 
is  given 

iZ/(aj  ,  y  )  =  -  (4/.V/)f  (cu)  y  ^  Y k sink ^sin  k sink ^  (10) 

m  n 

with  k^^  mil/ L^,  k^  =  nn/ Ly, 

(a^^,/i^f-(k^/kjy^2(k^/k^y(kjk,yf-(kjk^y)- 
=  [£'xy/(l  -  v^)Pp/i coy  =  CO 

c,  =  [£/(l-vypy‘^" 

r,  =  {I/hy'^ 

M  =  p  ptiL^Ly 

In  the  above,  pp  is  the  plating  mass  density,  I  its 
cross-sectional  moment  of  inertia  h  the  plate  thickness,  (CpT^)^ 
and  (CpT g)y  are  the  products  of  the  compressional  speeds  in  the 
plating  material  and  the  cross-sectional  radii  of  gyration,  along 
the  x,y  directions  respectively,  and  kf^  and  kf^  are  the 
associated  flexural  wavenumbers. 


For  the  distributed  pressure  p{uo-,  x  ^  ,y  Eq.  10  serves  as  an 
influence  function  and  the  displacement  becomes 

!l/(aj;.Y,y)  =  -  (4/M)  ^  ^  /l^„r„„sin(/c^x)sin(/c„y) 

m  n 

with 


(12) 


(13) 
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Having  detairmined  the  plate  displacement,  the  corresponding 
velocity  and  acceleration  may  be  computed  from 

v(<jo;x  ,y)  =  -iojw(Lfj\x  ,y)  (14a) 

and 

d(uo;x  ,y)  =  -<xi^w(<x>\x  .y)  (145) 

The  associated  flexural  stresses  including  those  evaluated  at  the 
outer  most  plating  fibers  are 

a,  =  [E^(h/2)/(l-v^)][o^uu/3l^vd^w/dy^]  (15a) 


a=[E..(h/2)/il-v"^)][a^w/dy^  +  w~w/dl] 


(155) 


Specifically,  consider  the  harmonic  pressure  component  of  an 
acoustic  plane  wave  of  magnitude  Pi  incident  on  our  plate  from  a 
simple  farfield  acoustic  source  at 

Pi(uo;x',y',0)  =  exp  [-i/csm0;(x' cos  <});  +  y 'sin  <!),)]  (16a) 


with 


P^  =  pQexp(ikP)/4nP  (165) 

where  R  is  range,  and  Q  the  strength  (volume  acceleration)  of  the 


source. 


The  wave  impinges  on  the  plate  (t=0)  and  reflects.  Using  the 
high  frequency  asymptotic  solution  /cly»l)  for  analytical 

simplicity,  this  reflected  wave  causes  a  pressure  doubling. 
Therefore,  substituting 

p(cj;x',y')  =  2pj(cu:x',y',o)  (17) 

in  Eq.  13, 


A^^  =  2P,  exp[-j(K^x'  +  Kyy')]sin(A:,„x')sin(/c„y')dx'dy' 


=  2PJJ^ 


(18) 


with 
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=  [A:^exp(-iK^I^/2)/(A:^-  K^)][exp(;K,^I^/2)-  cos(mn)exp(-iK_^Z_^/2)] 


/y  =  [fc„exp(-JKyi:y/2)/(;c^-  Ky)][exp(/KyZ.y/2)-  cos(mn)exp(-/Kyiy/2)] 


=  /c  sin  0,  cos  (j), 
K  =  ^ sin  0,  sin  (|), 


(19) 


The  impulse  response  corresponding  to  any  of  the  above  response 
functions  is  given  by  the  transform 


G(t-,x,y)  =  n~  Rq\  G((x>;  x  ,y)exp(^- L(X)t)ciuo 


(20) 


Using  Eq.  20  as  the  Green's  function  for  an  incident  wave  of 
arbitrary  waveform  p,(0  the  general  time  domain  solutions  for  our 
response  functions  may  be  computed  from  either  the  convolution 


t 

f(.t-,x,y)=j  p^(t')G(t-t';x,y)df 


(21a) 


or  the  inverse  transform 


f  (t-,  X  .y)  =  n'  Re\  P;(a))G(cu:x,y)exp(-iooOdoo 


with  the  frequency  spectrum  of  the  incident  wave  given  by 


(216) 


Pi(^)=  j  Pi{t)exp(i(Jot')dt  (22) 

0 

For  example,  with  a  classic  unit  amplitude  N-shaped  wave  of 
duration 

p,(0  =  ( 1 -2f/T  J  0<(<x„  (23a) 

the  spectrum  is 

p,(aj)  =  (tcu)'*[2exp(ia)Tp/2)sin(ujTg/2)/((juT„/2)-(l-*-exp(iUj'Cj,)] 

(236) 

although  there  are  no  such  restrictions  on  p,(0- 
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3.1.2 


Sound  Radiated  By  Mechanically  Driven  Plate 


Let  us  now  consider  the  sound  radiated  by  the  above  plate  when  it 
is  driven  mechanically,  as  opposed  to  acoustically  (Fig.  3b) . 

Eq.  10  still  holds  but  here  F (uj)  represents  the  spectrum  of  the 
direct,  or  mechanical  drive,  at  The  far  field  acoustic 

pressure  radiated  to  (/? ,  Gj ,  (p,)  by  this  vibration  field,  is^® 

=  -(oo’p/2ji)lv  (24) 

where  V',  the  double  spatial  Fourier  transform  of  the 
displacement  field,  is  defined  by 

W''((ju;k,.x.^)=  f  j  i2;(cu:.Y,y)exp[-(K,x  +  iCyy)]dxdy  (25) 

h  0 


In  Eq.  24  R  is  again  range,  K^^.Ky  are  as  defined  previously  in 
Eq.  19  and  as  before  for  consistency  (as  well  as  simplicity)  we 
have  again  invoked  the  high  frequency  approximation 
{kL^,  fciy»l),  here  by  embedding  the  plate  in  a  rigid  baffle  (The 
Kirchhoff  approximation) .  Combining  Eqs.  24  and  25 


p,((ju;i?,e,,<Pi)  =  -(ou^p/2ii)[exp(tA:;?)//?] 


1  y  i  y 

//  u)  ( (JO  ; -V ,  y  )  exp  [ -  j  ( ^  y  y  ) ] 


(26) 


or,  using  Eq.  10 

P,((J0:/?,  9,,  <j)i)/F((u;x,,y  J  =  (4/A'f)((u^p/2n)[exp(jfc/?)//?] 

iy 

yy^mnf  f  sin(A:^x')sin(/c„y')exp[-j(K^x'  +  Kyy')]rfx'dy'sin(fc,„xJsin(/c„y,) 
"  0  0 

(27) 

We  now  compare  Eq.  27  with  the  expression  previously  developed 
for  the  structural  acceleration  in  response  to  the  acoustic 
source  (Eq.  14b  using  Eq.  12  with  defined  by  Eq.  18  and  Ft  by 

Eq.  16b) 
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a  (  cu :  X  ,  y  )  /  Q  ( CAJ ;  ,  9 , ,  (t)  j  =  (  4  / .  W  )  ( .ju "  p  /  2  n )  [  e  X  p  ( i  A: .? )  /  ,?  ] 


V  V*  : 


I  sin(4^.x') 


sin(^'„y')exp[-^(K^x'+  K^y ')]d;c'dy'sin(,l-,^.Y)sin(4-^y)  (28) 

We  see  that  the  point  harmonic  acceleration  of  the  plate  at  (x,y) 


to  the  acoustic  source  Q  at  (.?.9..6.)  is  reciprocal  to  the 
pressure  radiated  to  by  the  plate  harmonically  vibrating 

in  response  to  the  mechanical  drive  F  at  (a'„  =  x .  y .  =  y ),  that  is, 

p,(aj;/?,e;,(j)^)/f(co;x,y)  =  d((X):x,y)/Q((jO,i?.9,.(l)J  (29) 


The  harmonic  point  displacements  and  velocities  mav  also  be 
determined  reciprocally  from  Eqs.  14.  Also,  in  light  of  Eg.  20 
reciprocity  holds  in  the  time  domain  for  impulses,  or  more 
generally  as  long  as  the  force  and  the  acoustic  source  exhibit 
the  same  spectra. 


Using  superposition,  reciprocal  relationships  may  also  be 
developed  for  distributed  rather  than  point  response  functions, 
for  example  the  average  acoustically  excited  plate  response  is 
reciprocal  to  the  radiation  from  a  uniformly  driven  plate. 
Reciprocal  relationships  also  exist  for  response  functions 
involving  gradients,  for  example  stress  (Eqs.  15).  However,  here 
reciprocity  requires  a  mechanical  drive  with  a  higher  order 
singularity,  e.g.,  a  doublet  for  a  moment,  etc.  This  is  deemed 
impractical . 


3.1.3  Numerical  Results 

In  this  section  we  present  numerical  results.  The  geometric, 
temporal  and  physical  constants  are  chosen  for  illustrative 
purposes  but  are  representative.  Unfortunately,  even  with  our 
elementary  model  there  is  a  multitude  of  parameters.  They  are 
nondimensionalized  as  follows  for  the  purpose  of  generating 
universal  response  curves 

y=Ly/Lx  =  ratio  of  plate  lateral  dimensions 
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r,/L,  =  Ch./-il2L,^  =  slenderness  ratio 
r\  =  material  loss  factor 
V  =  Poisson's  ratio 

(Cprg)^/(Cprg)y  =  ratio  of  plating  stiffnesses  along  x  and  y 


axes 

a  =  sin'‘(l/M)  =  9i  =  Mach  angle,  or  polar  angle  of  the 
incident  acoustic  wave  or,  reciprocally,  of  the  radiated 


pressure . 

M  =  c/i/  =  effective  Mach  number 

(}),  =  azimuthal  angle  of  the  incident  acoustic  wave  or 

reciprocally,  of  the  radiated  pressure 

Cp/c  -  ratio  of  plating  sound  speed  to  that  of  air 

=  cx^/L^  =  duration  of  sonic  boom  normalized  by  acoustic 
propagation  time  over  the  plate  span. 

(x / Ly ,y / Ly)  =  normalized  coordinates  of  observer  location  on 
plate  with  acoustic  excitation  or  reciprocally,  the 
mechanical  drive  location  on  the  radiating  plate 
kLx  =  <->^Lx/c  =  plate  lateral  dimension  along  x  axis 
normalized  to  acoustic  wavenumber  at  frequency  co  in  air 
=  (iuo/Cprg)Y^L^  =  plate  lateral  dimension  along 
X  axis  normalized  to  flexural  wavenumber  at  frequency  <x>  in 
plate 

x=t/Xg  =  time  normalized  to  pulse  duration 
w  =  w/Wst  =  plate  displacement  normalized  to  the  static 
displacement  of  (square)  plate  in  response  to  uniform 
pressure  (lo,,  =  0.049(1, //i)"(  1  - 

a  =  a/0st  =  flexural  stress  normalized  to  the  static 
stress  at  (square)  plate  center  in  response  to  the  uniform 
pressure  Pi  (0^,  =  0.22(Z,y//i)^(  1  +  v)/^ i) 

V  =  vx „/ w St  =  nondimensionalized  plate  velocity 
a  =  ax'lwst  =  nondimensionalized  plate  acceleration 


For  our  illustrative  example  we  take  Xg/L-lAxlO  ^  and 

(3  =  3.53x10"^  which,  for  example  models  a  perfect  N-wave  of  x^  =  0.2 
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sec  duration  incident  on  a  square  uniform  plate  2 . Sen  thick  with 
lateral  dimension  2.5m.  To  approximate  solid  wallboard 
construction  we  further  assume  c^/c  =  4.87 , v =  .25  and  somewhat 
arbitrarily  r|=.l.  For  perspective,  with  these  parameters 
11  "^0  =  3.08  and  therefore  the  ratio  of  the  period  of  the 
fundamental  natural  frequency  of  the  plate  to  the  pulse  duration 
is  7  11 /Tg  =  2n/coTo  ~  2.  To  determine  response  time  histories  the 
associated  frequency  spectra  have  been  computed  up  to  cjUTg  =  50  at 
an  increment  of  Aoutg  =  0.0  1. 

In  Fig.  4a  we  plot  the  computed  spectra  for  the  boom  itself  and 
the  acoustic  impulse  response  of  the  plate  displacement  at  the 
center.  Also  shown  as  a  reference,  is  the  spectral  displacement 
of  a  simple  distributed  oscillator  (one  degree  of  freedom)  when 
the  mass  is  driven  directly  by  the  harmonic  pressure  with 
magnitude  Pi. 

The  oscillator  frequency  is  equal  to  the  fundamental  plate 
frequency  and  its  response  has  similarly  been  normalized  to  its 
static  value.  The  low  frequency  plate  displacement  is  6  dB  above 
that  for  the  oscillator  due  to  the  assumed  pressure  doubling  with 
the  acoustic  excitation.  The  corresponding  acoustic  impulse 
spectra  for  the  plate  velocity,  acceleration  and  flexural  stress, 
all  evaluated  at  the  center,  have  been  plotted  in  Fig.  4b.  As 
with  the  displacement  the  (eventual)  6  dB  low  frequency  value  for 
the  stress  is  attributable  to  pressure  doubling.  At  the  higher 
frequencies,  and  ignoring  oscillations,  the  acceleration  is 
essentially  frequency  invariant. 

Of  particular  interest  are  the  velocity  and  stress  spectra.  As 
discussed  earlier  in  Section  2.1,  Hunt^  has  shown  that  for 
resonant  rods  and  low  order  resonant  beams  and  plates,  the  ratio 
of  peak  stress  to  peak  velocity  is  of  simple  form 

0  pk^ijpkP  pCp^K  (30) 
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Fig.  4  Sonic  Boom  and  Plate  Response  Spectra:  (a)  Sonic  boom, 

plate  and  oscillator  displacement;  (b)  Flexural  stress 
velocity  and  acceleration. 
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where  and  Cp  refer  to  the  density  and  sound  speed  in  the 
structural  element.  For  a  simply  supported  plate 

K  =  V3(  1  +  av)/4(  1  +  Qt^)  (31) 

With  our  parameters  (a=l,v=.25)  K=1.08  or  0.7 dB.  This  result 
is  explored  in  Fig.  5  where  the  spectra  of  the  stress  and  the 
ratio  of  stress  to  velocity  have  been  replotted,  the  latter 
normalized  as  follows 

5/i?PpCp  =  [4.49-j  12(rj/i^)(l  v)(Cp/c)/,3]a/y~2.98a/y  (32) 

or 

20  log  |a/iJppCp  |=9.5  +  201ogo-20Iogy 

Also  plotted  in  Fig.  5  is  Eq.  31.  As  expected,  the  exact 
solution  is  well  predicted  by  Eq.  31  at  the  lower  natural 
frequencies  of  our  plate  where  the  stress  peaks,  although  the 
comparison  may  be  poor  elsewhere. 

Computed  time  histories  are  shown  in  Figs.  6  and  7a  with  that  of 
the  boom  as  a  reference.  The  long  time  oscillations  (ringing)  by 
and  large  correspond  to  the  fundamental  plate  frequency.  Hunt's 
relationship  between  stress  and  velocity  is  once  again  verified 
in  the  time  domain  where  we  see  that  the  ratio  of 

=  2.476.3  =  0.38  or  -8.4  dB  and  therefore  Op^/fp*.ppCp~  1 . 13  or  1 
dB  (vs  0.7  dB) .  Finally,  in  Fig.  7b  we  show  stress  and  velocity 
and  time  histories  for  an  orthotropic  plate  and  oblique 
incidence/radiation.  Here  we  compute  Op^^/fp^PpCp-O.ZS. 

3 . 2  THE  DIFFRACTED  AIRBORNE  PATH 

In  practice  structural  elements  such  as  individual  windows  and 
walls,  may  often  be  shadowed  from  the  boom  and  also  baffled  by 
surrounding  structure.  Structural-acoustic  reciprocity  allows 
for  both  of  these  phenomena.  This  is  illustrated  below  for  the 
example  shown  in  Fig.  8  using  a  coupled  finite  element-boundary 
element  (FE/BE)  numerical  approach  to  the  solution  of  the 
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Fig.  6  Time  Histories  for 
Oscillator  Displace 


Fig. 


T=t/To 


7a  Time  Histories  for  Sonic  Boom  Pressure  and  Plate  Velocity 
and  Stress.  ^ 

(e.=(}).=0,  (r  )  /(r  )  =1:  (a/vp  c  )  =2.98  a/v~1.15 

max 
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Fig.  7b  Time  Histories  for  Sonic  Boom  Pressure  and  Plate  (center) 
Velocity  and  Stress. 

(»i=<l>i=45°,  (r  )  /(r  )  =4:  (o/vp  c  )  -2 . 98o/v~0 . 75 

^  y  ^  X  P  P  max 
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Helmholtz  integral  equation.  The  example  also  illustrates  at  the 
currently  available  power  of  state-of-the-art  general  purpose 
FE/BE  computer  codes  to  model  complex  structural  acoustic 
interactions . 

3.2.1  The  Helmholtz  Integral  Equation 

The  dynamic  behavior  of  a  structure  embedded  in  an  acoustic 
medium  can  be  described  by  the  matrix  equation!^ 

Zll  =  F -CAP  (33) 

where  Z  is  the  structure’s  mechanical  impedance,  u  is  the 
velocity  vector,  F  is  the  mechanical  force  acting  on  the  body,  G 
and  A  are  the  geometry  matrices  of  the  surface  normals  and  areas 
defining  the  fluid-structure  interface,  and  P  is  the  fluid 
pressure  acting  on  the  body.  The  structure's  mechanical 
impedance  Z  is  defined  as  a  function  of  frequency  as 

Z  =  F /u.  =  (i(x>M  +  D K /iuo)  (34) 

where  M,D,  and  K  are  the  mass,  damping  and  stiffness  matrices  of 
the  structure. 


The  fluid  behavior  is  described  by  the  Helmholtz  integral 
equation 


aZ)(r) 

PM---  ~  dS- 
—  dn 


q{x)D{r)cLS 


I p{x')/2-  Pi,  X  on  S 
\p(x')-p;,  X  on  E 


(35) 


where  the  vectors  x,  r,  x',  and  n  along  with  the  geometry  are 
described  in  Fig.  8,  Z?(^)  is  the  Green's  function,  and  q{x)  is 
the  normal  pressure  gradient.  These  variables  are  defined  in  the 
equations : 
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D{r')  =  e  ‘*-/4jtr , 


(36) 


g(x)  =  3p(x)/3n  =  -iiopd„(x) 


(37) 


d  D(r)/ dn^ - (Lk+  l/r)cos|3 

-  4itr 


(38) 


The  quantity  d„  in  Eq.  37  is  the  normal  velocity  on  the  surface, 

S,  while  the  far-field  angle  [3  relative  to  the  surface  normal  is 
defined  in  Fig.  8.  Substitution  of  Eqs.  36,  37,  and  38  into  the 
form  of  Eq.  35  valid  for  points  on  S  (taken  to  be  the  surface  of 
the  structure)  yields  the  integral  equation 

-p(x')-  I  p(x)^ - \  ik  + -  ]cosRdS  =  iodo  )  u„(x)- - dS  +  p,  (39) 


J  p(x)  — |^iA:  +  -Jcos(3dS  =  iaip  J  u„(x)  — dS  +  P;  (39) 

s  *  —  s  — 

Numerical  integration  of  the  surface  integrals  allows  Eq.  39  to 
be  cast  into  a  discrete  matrix  equation  of  the  form 

£f  =  Cu„  +  P..  (40) 

The  normal  velocities  on  the  structure  are  obtained  from  the  G 
matrix  by  the  equation 


U„  =  GT  LL. 


Eqs.  33,40,  and  41  can  be  combined  to  give 

HP  =  Q  +  Fi, 


(41) 


(42) 


where  H  =  E  +  CG^ Z'^GA  and  Q  =  GG^Z~^F. 


Eq.  42  is  the  coupled  matrix  equation  for  the  structure  and  the 
fluid  valid  on  the  surface  of  the  structure  which  is  solved  to 
give  the  pressure  acting  on  the  surface.  The  normal  velocity 
is  then  obtained  from  Eqs.  33  and  41.  Substituting  these 
quantities  into  Eq.  35  provides  the  solution  for  the  far— field 
pressure . 

Lastly,  we  note  that  the  Helmholtz  integral  equation  is  singular 
at  a  discrete  set  of  "forbidden"  frequencies.  These  frequencies 
correspond  to  the  acoustic  resonances  of  the  interior  cavity.  To 
alleviate  this  problem,  NASHUA  introduces  "CHIEF  points" 
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(combined  Helmholtz  Integral  Equation  Formulation^®)  to  the 
interior  volume.  These  points,  together  with  the  interior 
Helmholtz  integral  equation,  are  used  to  generate  a  number  of 
additional  equations  that  impose  needed  constraints  to  the 
system. 

3.2.2  Illustrative  Example 

Our  example  is  of  a  cubic  structure  resting  on  a  rigid  ground 
plane  (Fig.  9a).  All  walls  of  the  structure  save  one  are  assumed 
immovable.  Centrally  embedded  in  the  otherwise  rigid  remaining 
is  a  square  "window"  modelled  as  an  isotropic  simply 
supported  plate.  The  window  side  length  (1^,)  is  one  third  that 
of  the  cube  (i)  and  its  slenderness  ratio  is  =  4. 8x  10“^.  it 

is  made  of  glass  (Cp/c=  15.5,  v  =  0.3  and  somewhat  arbitrarily, 
ri=10  *)  .  For  computational  convenience  the  interior  air  space 
is  ignored.  The  finite  element  model  of  the  structure  is  created 
using  the  finite  element  program  COSMIC  NASTRAN.  The  formulation 
of  the  Helmholtz  boundary  elements  and  its  interface  with 
NASTRAN  is  performed  using  the  NASHUA  code^'^. 

To  numerically  demonstrate  reciprocity  the  window  will  be  driven 
acoustically  and  mechanically  and  the  appropriate  results 
compared  in  the  spectral  domain  for  nondimensional  frequencies 
0.1<kL<1.2.  A  computer  plot  of  our  model  is  shown  in  Fig. 

9b.  Because  the  structure  and  the  applied  loads  are  symmetric 
only  half  the  structure  is  modeled.  Highlighted  in  the  plot  is 
the  plane  of  symmetry,  the  ground  plane,  and  the  window.  Each 
element  measures  roughly  1/8  of  an  acoustic  wavelength  at  kL=  12. 
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Fig  9a.  Elastic  Window  Embedded  In  An  Otherwise 
Rigid  Walled  Cubic  Structure  Resting  On  A 
Rigid  Ground  Plane 


L/3 


3.2.3  Numerical  Results 


First,  for  our  reciprocal  calculation  we  apply  a  concentrated 
normal  force  to  the  center  of  the  window  and  compute  the  radiated 
pressure  at  a  number  of  receiver  locations  assumed  to  be  in  the 
far  field.  Specifically,  the  pressure  is  calculated  at 
(i?,0,0).  (;?  ,11/2,  Ji/2)  and  (/?,Jt/2,0)  corresponding  to 

receivers  in  front  of,  in  back  of,  above,  and  to  the  side  of  the 
window  and  where  is  a  reference  far  field  range.  Results 

are  plotted  in  Fig.  10.  As  a  reference  with  L  =  3L^  =  2.2Sm  and 
assuming  c  =  340m,/s,  kL=  1.0  corresponds  to  24.  Hz.  The  radiated 
pressure  has  been  no2rmalized  to  that  radiated  in  the  absence  of 
the  rigid  structure  but  in  the  presence  of  the  rigid  ground 
plane. 

At  small  values  of  kL  the  driven  window  behaves  as  a  monopole 
source  radiating  omnidirectionally.  At  the  receiver  located 
above  the  window  destructive  interference  from  ground  reflections 
is  evident  beginning  at  fcZ,  =  it  where  the  path  difference  between 
the  direct  and  the  ground  reflected  waves  first  measures  one-half 
acoustic  wavelength  in  air.  The  baffling  effect  of  the  rigid 
walls  surrounding  the  window  accounts  for  up  to  a  10  dB 
difference  in  front  to  back  radiation  levels.  Of  particular 
interest  is  the  fact  that  the  maximum  radiated  levels  actually 
exceed  the  6  dB  increase  associated  with  an  acoustic  source  in  a 
baffle  of  infinite  extent. 

For  our  direct  calculation  we  now  compute  the  velocity  at  the 
center  of  the  window  for  an  incident  plane  wave.  Corresponding 
to  the  earlier  receiver  locations,  the  orientation  of  our  plane 
wave  is  defined  by  a  far  field  acoustic  source  located  at 
spherical  coordinates  (/?0,0),  {R,n,0)  and  (^ ,  ti/2  ,  k/2),  the  fourth 

location  omitted  for  computational  convenience. 
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on  Window  Surface  Pressure 
ident  Pressure  Level. 


Direct,  and  reciprocal  calculations  are  compared  in  Figs.  11a, b 
and  c,  the  latter  renormalized  to  demonstrate  reciprocity  as  per 
Eq.  7  (or  29) .  Differences,  which  may  be  attributed  to  the 
discretization  process  and  the  accuracy  of  our  solution  to  the 
simultaneous  equations  given  by  Eq.  42,  are  generally 
imperceptible.  To  highlight  the  diffraction  effect  itself,  we 
plot  in  Fig.  12  our  direct  calculation  of  the  (normalized)  window 
center  velocity  for  the  ground  level  source  placements 
corresponding  to  marginally  supersonic  flyovers,  i.e.,  M~1 
oriented  such  that  the  window  is  illuminated,  shadowed  and  grazed 
by  the  incident  boom. 

Finally,  we  note  that  the  above  coupled  FE/BE  approach  can  have 
easily  been  used  to  analyze  a  more  complex  geometry  and 
structure.  For  example,  the  walls  need  not  have  been  rigid  nor 
in  the  shape  of  a  cube.  The  window  could  also  have  been  of  any 
shape  and  of  any  construction  allowing  linear  viscoelastic  or 
structural  modelling.  This  could  include  provisions  for  studded, 
composite,  or  brick  and  mortar  walls  and  brie  and  brae 
attachments . 

3.3  THE  SEISMIC  PATH 

A  sonic  boom  induced  acoustic  field  impinges  on  the  ground  plane 
as  well  as  on  man-made  structures.  Generally,  due  to  the  large 
impedance  mismatch,  this  field  is  almost  completely  reflected 
from  the  plane  with  little  energy  coupled  into  the  ground.  The 
exception  is  at  "coincidence”  where  the  wavenumber  of  the 
airborne  acoustic  wave  projected  on  the  plane  matches  that  of  a 
characteristic  wave  propagating  in  the  ground.  In  particular, 
coupling  may  be  strong  at  coincidence  with  Rayleigh  surface  waves 
which  propagate  at  close  to  the  ground  shear  speed  and  remain 
trapped  near  the  surface.  As  with  the  airborne  pressure  field 
this  coupled  ground-borne  energy  may  reradiate,  excite  structures 
and  potentially  cause  unacceptable  noise  levels  and/or  damage. 
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Transfer  Function,  dB 


Fig.  11c  Transfer  Function  For  Window  Response  to  Acoustic  Wave 
Computed  Directly  and  Reciprocally:  Wave  Incident  From 
Above  (R,Tr/2,Tr/2)  . 
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In  practice  however  this  path  appears  to  be  relatively  weak.^^ 
Nevertheless,  it  also  satisfies  reciprocity  as  is  demonstrated 
below,  assuming  the  incident  field  to  be  acoustic  and  the  ground 
a  linear  viscoelastic  medium. 


For  our  illustrative  purposes,  we  consider  the  situation  pictured 
in  Fig.  13a.  As  before,  a  distant  acoustic  source  of  volume 
velocity  Q,  or  volume  acceleration  Q,  generates  an  effectively 
plane  acoustic  wave  in  the  far  field.  Letting  (j),  =  0, 

P;(uj  :  X  ,  z)  =  P,(cjo  :x ,  2)exp{iA:[(.r  -  Xj)sin  0  +  (z  -  z^)cos9]}  (45) 

with 


P,(a):x,z)  =  -iujp(2(a3;Xj,2j/4n/?j(|x-xJ,  |  z-zj) 
Evaluating  Eq.  45  on  the  ground  plane 

p,(a):x,0)-=^  =  :x,  0)exp{f/c[(x  -  Xj)sin  9  -  z^cos9]>  (46) 

The  reflected  (plane)  wave  is  given  by^O 

p(a):x,z)/Pi(oo:x,O)  =  Yexp{t7c[(x-Xj)sin0-(z-Zj,)cos0]>  (47) 


with 


Y  =  (Z(fc^)cos0-  l)/(Z(fc^)cos9+  1) 


(48) 


Z(/cJ  =  Z(/cJ/pc  =  (p^c,/pc){[l-2(fc,/^',)"]"-4(fc,/fc,)(A;,/fcj2 
[  1 -(/c,//cj"]‘'"[  l-(/c,//c, )"]''">/[  l-(fcy A:, )"]‘^" 


(49) 


where 


k  ^  =  ksiu  Q 


k^  =  uo/ c, 


k^  =  (jo/c. 

In  Eq.  49  Pg  is  the  ground  density  and  and  are  the 
dilatational  and  shear  ground  speeds.  The  total  (normalized) 
pressure  is 

Pf(co:x,z)/7,(a):x,0)  =  [F,(a):x,z)  +  p(a):x,z)]/7,((ju;x,0)  (50) 

which  evaluated  at  z  =  0 
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FIs  13a.  Direct  Model  of  insonified  Elastic  Ground  Plane 


Covering  Elastic  Half  Space 


0  Q(a);x  ,z  ) 


\  \ 


9 


v{a);x,0) 


Ground  Plane 


Elastic  Half  Space 


Fig  13b.  Reciprocal  Model  of  Sound  Radiated  by 
Mechanically  Driven  Ground  Plane 
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{2Z(/l\)cos  0/[  Z(a:^)  cos9  +  1  ]}  exp  {i/c[(x  -  x  j)  sin  9  +  x^cos  0]}  (51a) 

-2exp{ifc[(x  -  x„)sin  0  -  x^cos0]}  Zcos0»l  (515) 

The  associated  normal  particle  velocity  at  ground  level  is 

pc!j(co:x,0)/Z,((A):x,0)  =  p,((ju;x,0)//’j((X):x,0)Z  (52a) 

or 

p  cty(cu  ;  X  ,  0)  = -/oup  Q(a) :  x^ ,  X  j)[cos9/(Z(/'c^)  cos9  +  l)]e  W2nR^  (525) 

Now  consider  the  reciprocal  situation.  As  is  the  case  in  a  free 
field,  acoustic  reciprocity  allows  one  to  switch  source  and 
receiver 

p t((x> ;  X ,  z) / Q ((X)  ■,  X ^ ,  z =  p .(m  ;  X ^ ,  z  J / Q (<jo  ■,  X ,  z)  (53) 

In  addition,  for  the  reciprocal  equivalent  to  the  acoustically 
induced  seismic  motion,  consider  Fig.  13b.  A  compact  normal 
force  (F)  is  applied  to  the  ground  surface  at  (x,0,0).  The 
associated  axisymmetric  surface  velocity  may  be  expressed  in 
terms  of  its  radial  wavenumber  spectrum^! 

eo 

V{iX);r)=  j  v(ud\k,)J ^(k,r)k,dk,  (54) 

0 

with 

v(ix-,k,)  =  (F/2n)/[Z{k,)^z,(k,)]  (55) 

where  the  wavenumber  transformed  acoustic  impedance  is  given  by 

z,(k,)  =  pc/[l-(k//kf]''^ 

This  transformed  velocity  also  defines  the  far  field  radiated 
pressure, ^2 

ik  R 

p(co  ; ,  0^)  = -oop  iJ(ou  ; /c  sin  0^)e  V/?, 

.  iup 

= -iojp  F/[Z(/c  sin  0)  +  1  ]g  '"/2jt/?,  (56) 

Comparing  Eqs.  52b  and  56  we  arrive  at  the  reciprocal 
relationship 

?;(cu  ;x.O)/0(a):F,  =  /?,)  =  a(cju:x.O)/0(cA):/?,  =  /?,)  =  p(ou:/?,)/F(oo:x,0) 

(57) 
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Finally,  from  Eq.  49  we  note  that  since  p^Cd/pc»l,  Z»1 
ground  is  essentially  impenetrable  except  when 

(58) 

This  is  the  Rayleigh  wave  coincidence  condition.  In  contrast,  at 
dilatational  wave  coincidence,  A:^//Cd=l,  the  impedance  becomes 
infinite  and  the  ground  is  again  impervious. 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  analytical  results  presented  herein  are  promising  and  suggest 
that  further  work  be  directed  toward  empirical  validation  of  the 
proposed  techniques.  Specifically,  based  on  this  study  we 
conclude  the  following 

1.  Structural-acoustic  tests,  with  stationary  sources, 
provide  a  viable  technique  for  measuring  transfer  functions 
relating  sonic  boom  overpressures  in  the  vicinity  of  a  structure 
to  its  response,  either  in  terms  of  displacements,  velocities  or 
accelerations . 

2.  Such  transfer  functions  may  be  used  to  perform  damage 
assessments  in  terms  of  measurable  quantities  in  the  absence  of 
dedicated  flyovers. 

3.  Structural-acoustic  reciprocal  tests,  whereby  structures 
are  mechanically  driven  and  the  radiated  pressure  monitored  by  an 
array  of  microphones,  may  be  used  to  supplement  direct  testing 
and  increase  the  practicality  and  efficiency  of  performing  the 
required  measurements. 

4.  The  structural-acoustic  transfer  function  concept, 
direct  or  reciprocal,  is  modular  and  allows  for  easy  synthesis 
and  updating  with  related  work,  for  example  improved  codes  for 
predicting  ground  level  sonic  boom  spectra. 

5.  General  purpose  structural-acoustic  finite 
element/boundary  element  (FE/BE)  computer  codes,  introduce  a 
powerful  relatively  new  tool  for  predicting  the  response  of 
structures  to  sonic  boom  overpressures  potentially  accounting  for 
both  geometric  and  structural  complexities.  (For  structures 
modelled  as  linear  systems  the  predictions  inherently  satisfy  the 
appropriate  reciprocal  relationships.) 
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6.  Damage  predictions  using  reciprocal  measurements  may  be 
quantitative  with  levels  of  confidence  comparable  to  those  with 
current  techniques.  Predictions  may  also  be  qualitative, 
providing  a  ranking  of  vulnerability  among  different  structures 
and  environments  or  among  potential  damage  mechanisms  within  a 
given  structure.  In  either  case  such  measurements,  perhaps 
supplemented  by  FE/BE  analyses,  provide  a  valuable  input  to  those 
responsible  for  site  planning  and  development  in  environments 
vulnerable  to  supersonic  operations. 

Consequently,  we  recommend  that  the  proposed  techniques  be 
pursued  further.  Specifically  we  recommend  the  following: 

1.  Both  direct  and  reciprocal  structural-acoustic  testing 
techniques  for  measuring  transfer  functions  relating  sonic  boom 
overpressures  to  structural  response  be  empirically  validated  by 
comparison  with  equivalent  data  taken  at  a  series  of  dedicated 
sites/structures  with  actual  flyovers. 

2.  Flyover  data  should  also  be  used  to  validate  predictions 
from  corresponding  detailed  FE/BE  structural-acoustic  numerical 
models. 

3.  Upon  successful  completion  of  items  1  and  2,  procedures 
be  codified  to  obtain  required  structural-acoustic  transfer 
functions  and  integrate  them  into  an  overall  damage  assessment 
package  that  includes  current  estimates  of  sonic  boom 
overpressure  spectra  from  flight  characteristics  and  propagation 
factors,  and  correlations  of  structural  response  functions  to 
damage  for  conventional  and  unconventional  structures. 
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